A Systematic Search for Corotating Interaction Regions in 
Apparently Single Galactic Wolf-Rayet Stars. II. A Global View 

of the Wind Variability. 



A.-N. Chene 

Departamento de Astronomia, Casilla 160-C, Universidad de Concepcion, Chile 
Departamento de Fisica y Astronomia, Facultad de Ciencias, Universidad de Valparaiso, Av. Gran Bretana 

1111, Play a Ancha, Casilla 5030, Valparaiso, Chile 
Canadian Cemini Office, Herzberg Institute of Astrophysics, 5071, West Saanich Road, Victoria (BC), 

Canada V9E 2E7 
acheneSastro-udec . cl 
and 

N. St-Louis 

Departement de Physique and Centre de Recherche en Astrophysique du Quebec, Universite de Montreal, 
CP. 6128, Succ. Centre-Ville, Montreal, Quebec, Canada H3C 3J7 
stlouisOastro . umontreal . ca 



ABSTRACT 

This study is the second part of a survey searching for large-scale spectroscopic variability 
in apparently single Wolf-Rayet (WR) stars. In a previous paper (Paper I), we described and 
characterized the spectroscopic variability level of 25 WR stars observable from the northern 
hemisphere and found 3 new candidates presenting large-scale wind variability, potentially origi- 
nating from large-scale structures named Co- rotating Interaction Regions (CIRs). In this second 
paper, we discuss an additional 39 stars observable from the southern hemisphere. For each star 
in our sample, we obtained 4-5 high-resolution spectra with a signal-to-noise ratio of ~100 and 
determined its variability level using the approach described in Paper I. In total, 10 new stars 
are found to show large-scale spectral variability of which 7 present CIR-type changes (WR8, 
WR44, WR55, WR58, WR61, WR63, WRIOO). Of the remaining stars, 20 were found to 
show small-amplitude changes and 9 were found to show no spectral variability as far as can be 
concluded from the data in hand. 

Also, we discuss the spectroscopic variability level of all single galactic WR stars that are 
brighter than v ^ 12.5, and some WR stars with 12.5 < v < 13.5; i.e. all the stars presented in 
our two papers and 4 more stars for which spectra have already been published in the literature. 
We find that 23/68 stars (33.8%) present large-scale variability, but only 12/54 stars (^22.1%) are 
potentially of CIR-type. Also, we find 31/68 stars (45.6%) that only show small-scale variability, 
most likely due to clumping in the wind. Finally, no spectral variability is detected based on the 
data in hand for 14/68 (20.6%) stars. Interestingly, the variability with the highest amplitude 
also have the widest mean velocity dispersion. 

Subject headings: stars: variables: other — stars: winds, outflows — stars: Wolf-Rayet 
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1. Introduction 

Up to date, no reliable direct observations of 
Wolf-Rayet (WR) rotation rates have been ob- 
tained. In fact, since the spectrum of WR stars 
originates exclusively in their hot and dense wind, 
the surface cannot be seen and no "classical" 
photospheric line can be observed. This renders 
the direct determination of rotational velocities 
of WR stars using Doppler broadening of photo- 
spheric lines impossible. However, in some cases 
the rotation period can be deduced for single WR 
stars presenting periodic wind variability. When 
present in a WR wind, this type of variability 
could be due to the presence of Co-rotating Inter- 
action Regions (CIRs) generated by perturbations 
at the base of the wind which propagate through 
it and are carried around by rotation. Although 
the detected period is not necessarily directly the 
rotation period of the star, there is likely a link 
between both periods that can be untangled with 
a proper monitoring and modeling of spectral lines 
(e. g. the case of the B su pergiant star discussed 
in lLobel fc Blomme llioosl ). 

Only two WR stars were known to show CIR- 
type variability; i.e. WR6 and WR134 which 
have a unique photometric period also present in 
the line-profile variability ( Firmani et al. I llQSOt 



(jAntokhin et al. I Il995t iMarchenko et al. \ \l99§t) . 
WC9 stars were also excluded since many WC 
stars of this sub-type have be e n shown to be 



Morel et al.lll997t iMorel et al. Ill999[ ). A binary 
scenario has also been proposed to explain the 
periodic variations, but no clear evidence has 
ever been found for the presence of a compan- 
ion. In order to identif y more candidates, we 
(in St-Louis et al. 20091 hereafter referred to as 
Paper I) have obtained 4-5 spectra well spaced 
in time of 25 northern, galactic WR stars which 
are thought to be single and brighter than v ^ 
12.5. Those observations allowed us to establish 
a way to group the WR stars in terms of the 
strength of their spectral variability, i.e. stars 
that do not show variability based on the data 
in hand (NV), stars that show small-scale vari- 
ability (SSV) and stars that show large-scale vari- 
ability (LSV). Our CIR-type candidates are taken 
among the LSV group. However, we did not 
consider that all LSVs are good CIR-type vari- 
able candidates. Indeed, we found that WN8 
stars show variability that appears different in 
nature from those we are attempting to identify 
(see Paper I), and might have a different origin 



binar ies ( Hackwell et al. ] 119761: IWilliams et al. 



19781 ) and present spectral variations that, with- 
out proper monitoring, might look like CIR-type 
changes, but in reality be binary-related. These 
two classes of stars deserve a more thorough and 
dedicated study. Finally, follow-up observations 
of one candidate, WRl, have recently confirmed 
the presence of CIR-type variability in that star 



(IChene fc St-Louis 11201(11) 



In this paper, we complete the survey initiated 
in Paper I of a ll WR stars in the catalogue of 
van der Hucht I (|200lr ) brighter than v ~12.5, and 
we add stars down to a magnitude of ^ 13.5 for 
a total of 39 apparently single southern, galactic 
WR stars. We have not included any WN8 stars 
in this southern sample for the reason mentioned 
above. Also, some stars brighter than the cho- 
sen magnitude limit have been excluded from our 
list, i.e. WR25 for which an orbital pe riod has 
since been identified ( Gamen et al. 2006[ ). WR94 
which turned out to be a visual binary and WR 7 
and WR60, for which we unfortunately did not 
get enough spectra to investigate their variability. 
Our targets are listed in Table [T] where the stars' 
name, spectral type, RA, DEC, v magnitudes and 
terminal velocity, Voo (when available in the litera- 
ture), are provided. The last column indicates the 
stars' variability status (see below). 

Section [2] summarizes the observations. The 
search for variability is described in Section [3] and 
a discussion of the results for the new sample of 
stars and a summary of the spectroscopic variabil- 
ity level of all the Galactic WR stars in our global 
sample are presented in Section |4l 

2. Observations 

We carried out spectroscopic observations us- 
ing the 1.5m-telescope of the Cerro Tololo Inter- 
American Observatory (CTIO) during 14 nights 
between 2004 June 28 and July 12. As for the 
stars discussed in Paper I, we obtained 4-5 spec- 
tra per target with the long-slit spectrograph. We 
have reached a typical signal-to-noise ratio of 100 
over a spectral range of 5200-6000 A and with a 
resolution of 2.3 A (3 pix.). 

The data reduction was performed using the 
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standard IRAF software packages. First, a bias 
frame was subtracted, then each image was di- 
vided by a flat field. After extracting the spectra, 
the wavelength calibration was carried out using 
a He-Ar calibration lamp. As our spectra are not 
photometrically calibrated, the final step of our 
data reduction procedure was to rectify the spec- 
tra. To do so, we have fitted a low-order Legendre 
polynomial to wavelength regions with no strong 
spectral lines and divided our spectra by the fitted 
curve. 

2.1. Choice of lines 

The spectral interval chosen for this work is dif- 
ferent from that adopted for the stars presented 
in Paper I. Actually, in the previous work, the ini- 
tial wavelength interval ranged from 4500 to 6000 
A. However, due to problems with the wavelength 
calibration, the reddest half of the spectrum had 
unfortunately to be excluded and the useful spec- 
tral range was AA=4500-5200A. For WN stars, we 
were therefore forced to limit our study to the 
HeiiA4686 and HeiiA4860 lines, which are very 
similar. For WC stars, except for the very strong 
CiiiA4650 line, which is often blended with a few 
faint lines, only weak lines were available. Never- 
theless, the above-mentioned transitions were able 
to provide us with a diagnostic of the line-profile 
variability level in most cases. However, for this 
study, we prefered to choose a different wavelength 
interval in order to monitor lines that are poten- 
tially of greater interest (i.e. formed at different 
radii, more sensitive to density variability, etc.). 

In WN spectra, there are two relatively strong 
isolated lines which are interesting to study within 
the 800A-wide spectral interval available at the 
CTI0-1.5m telescope. These are HeiiA5411, of- 
ten among the strongest emission lines in the op- 
tical band, and HeiA5876, a flat-topped line only 
slightly perturbed by two NaD interstellar lines. 
These two emission lines of different ionization 
states do not have the same shape, which in- 
dicates that they are formed in different veloc- 
ity regimes of the win d , i.e. at different radii. 



Schulte-Ladbeck et al. I (|1995l ) present a plot of 
the ionization stratification of their model of the 
WN4 star WR6; in general Hei lines are found 
to form at higher velocity than Heii lines. For 
WC spectra, we chose CiiiA5696, a flat-topped line 
particularly sensitive to density variations, and 



CivAA5802/12, a doublet that is unfortunately 
blended in some cases with several other lines. 
For WN/WC stars, we looked at the HeiiA5411 
and ClvAA5802/12 lines. The two lines are also 
formed in different (althou gh partly overlapping) 
velocity regions of the wind. iDessart et al. ( 2000l ) 
for example, present line- formation regions for the 
WC8 star WR135; they find that the CiiiA5696 
line peaks at a distance about twice as far as the 
CivA5802/12 doublet. 

3. Variability Search 

To search for significant line-profile variability 
of the targets, the Temporal Variance Spectrum 
(TVS) of each dataset was calculated using the for- 
malism of Fullerton et al. (1996) and the quantity 



Sj(99%) 



(TVS), 



was obtained, where co 



_i(99%) 

is the reciprocal of the rms of the noise level in 
the continuum in a timeseries of N spectra. The 
value of Ej(99%) quantifies the level of variabil- 
ity at each wavelength: a spectrum that reaches a 
value of n varies with an amplitude n times higher 
than the variability measured in the continuum 
(which is assumed to be pure noise) with a confi- 
dence level of 99%. The spectrum of a given star is 
considered significantly variable at a given wave- 
length j if the value of I]j(99%) is significantly 
greater than 1. When a line is identified as sig- 
nificantly variable, it is possible to calculate the 
amplitude of its variability relative to its inten- 
sity. To do so, we have calculated in Paper I for 

1 /2 

ea ch wa velength j a modified TVS^ as defined 
in IChene (2007) and divided it by the line fiux 
{Sj — 1), where Sj is the weighted mean flux at 
wavelength J. This ratio is named the cr-spectrum. 
One should note that the calculation of a does not 
take into account instrumental variations due to 
the noise level. Thus, when the variation level of 
a given line is too close to the noise, which is the 
case for weaker lines, a is artificially high. That 
is why we manually set aj =0 when the variability 
at a given wavelength j is not clearly significant 
according to the value of I]j(99%) or when the 
line has a relative intensity lower than 1.3, which 
is the intensity limit of a spectral line from which 
we can investigate the variability for the present 
dataset. Using the S(99%) spectrum and the cr- 
spectrum, three main categories were defined in 
Paper I: stars showing no profile variability (NV : 
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E(99%) spectrum < 1), stars showing small-scale 
profile variability (SSV : cr < 5 %) and stars show- 
ing large-scale profile variability (LSV : a > 5%). 
One has to note that in the case of many WC stars, 
the global line intensity seems to vary from spec- 
trum to spectrum. These changes are specially 
big in the case of the broad, strong emission lines 
(with relative intensity greater than ^ 5), prob- 
ably due to the difficulty to fit the continuum in 
these regions. In these cases, we have normalized 
the line intensity in order to isolate the variability 
cause by bumps. 

For this study, we do not use the same spec- 
tral lines as in Paper I and one may wonder if this 
could affect our variability diagnostic. In fact, it 
has been shown in Paper I that the cr-value does 
not change significantly from line to line, as long 
as the noise level is not of the same order or higher 
than the amplitude of the variability. To demon- 
strate that this is also true for the lines we have 
chosen in this study, we re-observed two targets 
already presented in Paper I. In Figure [Ija), we 
present in the upper panel, a montage of residuals 
obtained by subtracting the mean from individ- 
ual spectra for WR115 (WN6). These observa- 
tions were obtained with the OSIS spectrograph 
at the Canada-France-Hawaii telescope (CFHT) 
during 6 nights in 2003 August (AA = 1.5 A (3 
pix.), AA = 4400-6100 A, S/N =150). Using the a- 
spectrum shown in the middle panel of Figure[lja) 
and the same criteria established in Paper I, we 
find that, as in Paper I, WR115 is designated 
as LSV. This is also the case for WC stars, as 
shown in Figure [ijb), where we present the mon- 
tage of residuals of WR 111 (WC5) obtained with 
the long-slit spectrograph of the Observatoire du 
mont Megantic (OMM) in 2005 and 2006 (AA = 
1.6 A (3 pix.), AA = 4500-6100 A, S/N =100). In- 
deed, this figure shows that, as in Paper I, WR 111 
is designated as NV (or on the verge of being clas- 
sified as SSV when looking at the CivAA5802/12 
line). Moreover, for the strongest lines, when vari- 
ability is detected, the diagnostic deduced from 
the cr-spectrum does not depend on the S/N of 
the spectra. 

Interestingly, Figure [T] also shows that even 
when the lines are formed at different radii in 
the stellar wind and in different physical condi- 
tions, the observed variability does not seem to 
change significantly in amplitude. Of course, only 



a better monitoring would allow us to determine if 
the spectral features causing the variability have 
different shapes depending on the emission line 
and/or if time delays exist between the variations 
observed in different lines. 




Fig. 1.— a) Top: Montage of the HeiiA5411 
(left) and IIeiA5876 (right) residuals (individual 
spectra - mean) for WR115 (WN6). For both 
cases, the scale factor of the ordinate is indicated 
in the top right-hand corner of the plot. HJD - 
IIJD(2000) is indicated on the y-axis. Second from 
top: mean spectrum. Second from bottom: cr- 
spectrum. Bottom: S (99%) spectrum, b) Same 
as a) for WRlll (WC5). 

Figures [2] to [21] present montages of residuals 
for all our targets. The graphs are organized by 
spectral type starting with WN stars from early to 
late types (Figures [2] to [121) followed by WN/WC 
stars (Figures [T^] and [T51) and finally WC stars 
from early to late types (Figures [T4| to [211) . In 
each of these figures, the cr-spectrum and the E 
(99%) spectrum are presented in the middle and 
the bottom panels respectively. In the following 
subsections, we will discuss in turn the variability 
for stars of different spectral types. Our conclu- 
sions concerning the spectroscopic variability level 
of each star is given in the last column of Table [1] 
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Table 1 

Our Sample of Single WR Stars in the Southern Hemisphere 



M 

IN 3,H1C 


opcc"Ci3ji xypc 


RA (2000) 


DEC (2000) 




't;oo(kms ^) 


V 3.1 13.D1111jy 


WR7a 


WN4/WC 


07 


20 


22.38 


-23 


43 


57.6 


(12.5) 






NV 


WR8 


WN7/WCE+? 


07 


44 


58.22 


-31 


54 


29.6 


10.48 


1505^ 




LSV 


WRIO 


WN5ha 


07 


59 


46.25 


-28 


44 


03.1 


11.08 


1375^ 


1100^ 


NV 


WR14 


WC7+? 


08 


54 


59.17 


-47 


35 


32.7 


9.42 


2095^ 


2900" 


ssv 


WR15 


WC6 


09 


13 


11.77 


-50 


06 


25.6 


11.72 


3600" 




NV 


WR17 


WC5 


10 


10 


31.92 


-60 


38 


42.4 


11.03 


2200" 




SSV 


WR18 


WN4 


10 


17 


02.28 


-57 


54 


46.9 


11.11 


1800^ 




NV 


WR21a 


WN6+0/a 


10 


25 


56.49 


-57 


48 


44.4 


12.80 






NV 


WR23 


WC6 


10 


41 


38.33 


-58 


46 


18.8 


9.67 


2240^ 


2900" 


SSV 


WR24 


WN6ha 


10 


43 


52.27 


-60 


07 


04.0 


6.49 


1985^ 


2160^ 


NV 


WR28 


WN6(h)+OB? 


10 


48 


58.68 


-59 


03 


37.5 


12.98 


1200^ 




SSV 


WR33 


WC5 


11 


00 


00.72 


-57 


48 


59.5 


12.35 


4500" 




NV 


WR44 


WN4+OB? 


11 


16 


57.86 


-59 


26 


24.0 


12.96 


1400^ 




LSV 


WR50 


WC7+OB 


13 


18 


01.07 


-62 


26 


04.5 


12.49 


3200" 




SSV 


WR52 


WC4 


13 


18 


28.00 


-58 


08 


13.6 


9.86 


2765^ 


3600" 


NV 


WR53 


WC8d 


13 


30 


53.26 


-62 


04 


51.8 


10.88 


1700" 




SSV 


WR54 


WN5 


13 


32 


43.79 


-65 


01 


27.9 


12.99 


1500^ 




SSV 


WR55 


WN7 


13 


33 


30.13 


-62 


19 


01.2 


10.87 


1200^ 




LSV 


WR57 


WC8 


13 


43 


16.37 


-67 


24 


04.9 


10.02 


1750^ 


2100" 


SSV 


WR58 


WN4/WCE 


13 


49 


04.52 


-65 


41 


56.0 


13.05 


1480^ 




LSV 


WR61 


WN5 


14 


13 


03.53 


-65 


26 


52.7 


12.41 


1400^ 




LSV 


WR63 


WN7+OB 


14 


50 


58.31 


-59 


51 


26.7 


12.83 


1700^ 




LSV 


WR67 


WN6+OB? 


15 


15 


32.61 


-59 


02 


30.8 


12.12 


1500^ 




SSV 


WR71 


WN6+OB? 


16 


03 


49.35 


-62 


41 


35.8 


10.23 


1200^ 




SSV 


WR75 


WN6 


16 


24 


26.23 


-51 


32 


06.1 


11.23 


2300^ 




SSV 


WR77 


WC8+OB 


16 


41 


19.12 


-48 


01 


59.6 


13.00 


2300" 




SSV 


WR78 


WN7h 


16 


52 


19.25 


-41 


51 


16.2 


6.61 


1335^ 


1385^ 


SSV 


WR79a 


WN9ha 


16 


54 


58.51 


-41 


09 


03.1 


5.29 






SSV 


WR79b 


WN9ha 


16 


55 


06.45 


-44 


59 


21.4 


8.32 






NV 


WR81 


WC9 


17 


02 


40.39 


-45 


59 


15.5 


12.71 


1300" 




LSV 


WR82 


WN7(h) 


17 


04 


04.61 


-45 


12 


15.0 


12.41 


1100^ 




SSV 


WR83 


WN5 


17 


10 


54.50 


-46 


36 


18.0 


12.79 






SSV 


WR85 


WN6h+OB? 


17 


14 


27.13 


-39 


45 


47.0 


10.60 


1400^ 




SSV 


WR87 


WN7h+OB 


17 


18 


52.89 


-38 


50 


04.5 


12.59 


1400^ 




SSV 


WR88 


WC9 


17 


18 


49.50 


-33 


57 


39.8 


13.25 


1500" 




LSV 


WR90 


WC7 


17 


19 


29.90 


-45 


38 


23.8 


7.45 


2100^ 


1800" 


SSV 


WR92 


WC9 


17 


25 


23.15 


-43 


29 


31.9 


10.43 


1130^ 


1300" 


LSV 


WRIOO 


WN7 


17 


42 


09.77 


-32 


33 


24.7 


13.44 


1600^ 




LSV 


WR103 


WC9d+? 


18 


01 


43.14 


-32 


42 


55.2 


8.86 


1060^ 


1300" 


SSV 



"^Spectral types are from Ivan der Hucht I (1200111 
HPrinia et aL I lll990f ) 
IHamann et aLI ll2006f) 
iTorres et aL I lll986fl 
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Fig. 2.— Same as FigureUJa) for a) WR 18 (WN4) 
and b) WR44 (WN4+0B?). 



Fig. 3.— Same as Figure for a) WRIO 

(WN5ha) and b) WR54 (WN5). 
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Fig. 4.— Same as FigureUfa) for a) WR61 (WN5) 
and b) WR83 (WN5). 
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Fig. 5. — Same as Figure [T^a) for a) WR21a 
(WN6+0/a) and b) WR24 (WN6ha). 



Fig. 6.— Same as Figure W[a) for a) WR28 
(WN6(h)+0B?) and b) WR67 (WN6+0B?). 




Fig. 7.— Same as Figure for a) WR 71 

(WN6+0B?) and b) WR75 (WN6). 
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Fig. 8.— Same as Figure [IJa) for a) WR85 
(WN6h+0B?) and b) WR55 (WN7). 
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Fig. 9.— Same as Figure [^a) for a) WR63 
(WN7+0B?) and b) WR78 (WN7h). 



6 




Fig. 10.— Same as Figure [I^a) for a) WR82 
(WN7(h)) and b) WR87 (WN7h+0B?). 



Fig. 13.— Same as Figure Ufa) for a) WR8 
(WN7/WCE+?) and b) WR58 (WN4/WCE), 
but for the HeilA5411 (left) and CivAA5802/12 
(right) line-profiles. 
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Fig. 11.— Same as Figure UKa) for a) WRIOO 
(WN7) and b) WR79a (WN9ha). 




Fig. 14.— Same as Figure UK b) for a) WR52 
(WC4) and b) WR17 (WC5). 



Fig. 12.— a) Same as Figure [Ha) for WR79b 
(WN9ha). b) Same as Figure HJa) for WR7a 
(WN4/WC), but for the HeiiA5411 (left) and 
CivAA5802/12 (right) line-profiles. 
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Fig. 15.— Same as Figure [UJb) for a) WR33 
(WC5) and b) WR15 (WC6). 
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Fig. 16.— Same as Figure [DJb) for a) WR23 Fig. 18.— Same as Figure [^b) for a) WR53 

(WC6) and b) WR 14 (WC7+?). (WC8d) and b) WR57 (WC8). 
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Fig. 17.— Same as Figure [^b) for a) WR50 
(WC7+0B?) and b) WR90 (WC7). 

3.1. WN stars 

A total of 5 WN stars are designated as NV. In 
the case of WRIO, WR21a, WR24 and WR79b 
this is simply because the emission lines that are 
present in the selected wavelength range are too 
weak (or simply absent) and no variability can 
be detected with the data quality in hand. How- 
ever, in the case of WR 18 (and to a lesser extent 
WRIO), in spite of the fact that the HeiiA5411 
line is quite strong, there is clearly no significant 
variability detected, excluding a few narrow peaks 
in the S-spectrum which barely exceed 1. 

All other WN stars are variable at either a large 
or small relative level. Interestingly, most of these 
stars, if not all, show significant variability in the 
absorption component of the HeiA5876 P Cygni 
profile. Moreover, the deeper the absorption, the 
higher the level of the E-spectrum. This feature 
is in fact the only variable part of the spectrum 
for WR28 and WR79a, thereby allowing us to 
classify them as SSV and not NV. Since it is not 
possible to clearly separate the absorption and the 
emission components of P Cygni profiles, and due 



Fig. 19.— Same as Figure [^b) for a) WR77 
(WC8+0B?) and b) WR81 (WC9). 




Fig. 20.— Same as Figure [^b) for a) WR88 
(WC9) and b) WR92 (WC9). 
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Fig. 21.— Same as Figure mb) for a) WR103 
(WC9d). 
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to blending with the CivAA5802/12 line, it is not 
possible to determine the exact depth of the ab- 
sorption component of the HeiA5876 line, and no 
relevant value of a can be computed. Therefore 
(7j=0 has been imposed at all wavelengths j cor- 
responding to an absorption in the mean profile. 
Note that WR 18 and WR 79b, both classified NV, 
do not show a clear P Cygni absorption in the Hei 
line. This is probably due to strong blending with 
the CIVAA5802/12 line in the case of WR18. The 
presence of stronger variability in the absorption 
components relative to the emission components 
of P Cygni profiles has already been mentioned in 
Paper I. We believe this is because the absorption 
part of the line profile comes from a much smaller 
volume of the wind than the emission part and 
therefore suffers less from cancelling effects from 
large-scale changes in density or ionisation arising 
in different parts of the wind. 

In addition to WR 28 and WR 79a, 9 other WN 
stars are classified as SSV. For WR54, WR67, 
WR71, WR75, WR82, WR83 and WR85, vari- 
ability has been detected at a level of 3-5 % of the 
line flux. The nature of the variability observed 
in these stars is similar to what was described in 
Paper I in the case of SSV stars. These changes 
are associated with clum ping in the WR win d, as 
observed for examp le bv lMoffat et al. ( 1988[ ) and 
Lepine et al.l ( 19961 ). The nature of the changes is 
less clear for WR78 and WR87. Indeed, in those 
cases, the I](99%)-spectrum shows only a few nar- 
row peaks which are not very significant and the o- 
spectrum has different amplitudes in the Hei and 
Hell lines. This is due to the low intensity of the 
emission lines that do not exceed 1.3 times the 
continuum. In the case of such a low-intensity 
hne, (Sj — 1) takes values close to and, when 
the variability has an amplitude as low as 2-3 % 
of the line intensity, the noise level constitutes a 
non-negligible fraction of the cr-spectrum. Conse- 
quently, cr reaches artificially high values. Note 
that for WR 71 and WR 75 the value of u reaches 
an amplitude of ^5 % of the line intensity which 
we adopted as the approximate limit between a 
SSV and a LSV classification. In this case, we de- 
cided to classify them SSV since the cr-spectrum is 
similar to what is observed for WRllO, designated 
as SSV in Paper I. 

When the line intensity is sufficiently high, 
large-scale variability levels are observed in WR44, 



WR55, WR61, WR63 and WRIOO. In all cases, 
we see a highly significant level of variability that 
reaches ~10 % of the line intensity with a a- 
spectrum that resembles the one presented in Pa- 
per I for WR134, a well-known CIR-type variable 



(iMorel et al.1ll999l) 



3.2. WN/WC stars 

In our sample, only three stars have a hybrid 
WN/WC spectral type. WR7a, is classified as 
NV because no significant variability is seen in 
the E(99%)-spectrum calculated for the moder- 
ately intense lines. WR8 and WR58 are classi- 
fied LSV. Interestingly, the spectral lines of WR 58 
have a globally variable intensity as can be seen 
from the very broad features in the residuals. This 
is probably due to variability in the continuum 
which translates as a variable dilution level of the 
spectral lines, since we have normalized the con- 
tinuum to one. Both WR8 and WR58 present a 
large-scale variability similar to what is expected 
for CIR-type variables. 

3.3. WC stars 

For WR15, WR33 and WR52, no significant 
variability is observed, even in very strong lines, 
and therefore, we classify them as NV. As for 
WR14, WR17, WR23, WR50, WR53, WR57, 
WR77, WR90 and WR103, the variability de- 
tected reaches a level of a=l-A % of the line inten- 
sity and therefore we classify them as SSV. Note 
that the high value of a in the Cm line of WR 17 is 
disregarded in view of its low intensity (see above). 
Also, without the normalization of the line inten- 
sity described above in section [3l WR 77 would 
show a much higher level of variability due to what 
seems to be real changes in the continuum. This 
is not the type of structure we are looking for and 
we prefer to keep it has SSV, but it would still be 
an interesting case to investigate. 

The only stars classified as LSVs are three 
WC9s out of four in our sample. Indeed, they 
reach a level of variability of more than 5 % of the 
line intensity. However, as described in Paper I, 
WC9 stars are possible dust producers although 
in this case they have not yet been identified as 
such. We will come back to this point in the next 
section. 
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4. Discussion and Conclusions 

This study completes our systematic search for 
Corotating Interaction Regions candidates among 
the brightest, presumably single, Galatic WR 
stars. In this section, we present a gfobal sum- 
mary of the results for the sample presented here, 
combined with those obtained for our northern 
sample and presented in Paper I. 

Following the method we have used in Paper I, 
we based our distinction between the different 
types of variability solely on the level of vari- 
ability. However, we could also use the velocity 
dispe rsion of the structures in the residuals. In- 
deed. iLepine fc Moffat I (|l999t ) define (^(cr^))LPv, 
i.e. the mean wavelet power spectrum of the 
residuals averaged over the spectral domain ^, 
to determine the level of variability. The maxi- 
mum of this power spectrum, ct^, gives the mean 
line-of-sight velocity dispersion of line profile vari- 
ability subpeaks. We c alculated (^((7f))T^ pv fol- 



lowing the formalism of Lepine fc Moffat ( 19991 ) 



for our complete sample of stars. In Figure [22l 
we present our results for WR18 (NV), WR50, 
WR71 and WRlll (SSV) and WRIOO, WR115 
and WR134 (LSV). The profiles for WRlll and 
WR134 are very simila r to those in Figure 4 of 
Leoine fc Moffat I (Il999h for the same stars, but 
with a spectral resolution more than ten times 
lower and a signal-to-noise of 100 instead of 200, 
we do not resolve as well structures with veloc- 
ity dispersion lower than 40 fcms"^. Moreover, 
the noise pattern can sometimes mimic true sig- 
nal at a few hundreds kms~^, and 5 spectra is 
not enough to average this artifact out from the 
mean wavelet power spectrum. To minimize these 
effects, we filter out any structures smaller than 
10 A;ms~^ and bigger than the terminal velocity 
Voo from the residuals. However, in the cases 
of low or very low level of variability, the mean 
wavelet power spectrum reaches the maximum at 
very large scales, which is meaningless when no 
significant variability is detected (e.g. WR18 for 
which we get ctj ~ 300 fcms ^). From Figure [22l 
we get CTj ~ 100 fcms^^ for the SSV cases and 
= 300 - 500 km s'^ for the LSV cases. To ver- 
ify if the distinction between the different types of 
variability is equivalent, regardless of the method 
used, we plot in Figure [23] the a^-value of the 
stars in our sample as a function of ((t(%)), i.e. 



the mean value of their cr-spectrum. Dotted lines 
mark the limit of cr > 5% adopted in this study, 
as well as ct^ > 250fcTOs~^, which seems to sep- 
arate SSV to LSV case s in both this study and 
Lepine fc Moffat I (|l999f ). Obviously, all the stars 
for which no good values of cr(%) could be com- 
puted due to a too low level of variability are ex- 
cluded from the plot. It is clear from this figure 
that stars with a high mean line-of-sight veloc- 
ity dispersion also have a high level of variabil- 
ity and that LSV are mostly located in the upper 
right quadrant of the plot. Some cases seem to 
disagree with the general trend. Indeed, WR90, 
WR 131 and WR 54 show variations with a level of 
variability significantly lower than 5%, but with a 
mean velocity dispersion higher than 250 fcms~^. 
But after inspection, we find that this is caused by 
the small number of spectra, since the stochastic 
appearance of clumps on the spectral lines com- 
bined with the noise pattern giveS the impression 
of shallow large-scale bumps on the line profile and 
these are clearly different from the structures we 
are looking for. As for WR63, high cr- values are 
obtained, even if the mean velocity dispersion of 
the variabilities is low. This is either due to the 
low intensity of the observed spectral line com- 
bine with a fairly high variability level (but still 
lower than 5%), which artificially boosts the a- 
value, or to the presence of real structures in the 
wind that actually have a small mean velocity dis- 
persion value. Because we cannot distinguish the 
two possibilities, we still keep WR63 in our list 
of LSVs, but we must keep in mind that this may 
not be the best CIR-type variable candidate. Fi- 
nally, we conclude that the determination of the 
type of variability is still more robust and easier 
to obtain when based on the level of variability, 
since this method is less dependent on the signal- 
to-noise ratio or on the spectral resolution than 
the method using the mean velocity dispersion of 
the variabilities. 

A compilation of the variability diagnostics as 
a function of the spectral type is presented in 
Figure [24l We add to this sample dat a from 
the literature f or the stars WR6 (Morel et al 

I995I) 



WR46 



19971). WR40 (lAntokhin et al. ■ 

(IVeen et al. 1 I2OO2I) and WR 138 fLepine et al 
19961) . 



Among WNE stars (including WN6 stars with- 
out hydrogen), 4 are classified NV, 10 SSV and 
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1 2 3 1 2 3 

log [cff (km s"')] log W( (km s"')] 

Fig. 22. — Mean wavelet power spectrum 
(ii(cr{))LPv for the spectral time series of seven 
stars in our sample. In the left panel are shown 
our results for WR18 (NV), WR50, WR71 and 
WRlll (SSV), which show a maximum at ~ 
100 km s ^. In the right panel are presented the 
stars WR 100, WR 115 and WR 134 (LSV), which 
show a maximum at ct^ 300 — 500 km s ^ . 
(i? (cTj)) LP V for WRlll was repeated in the right 
panel for better comparison. 



600 



500 r 




F = 

0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 
«7 (%)> 

Fig. 23. — Maximum of the mean wavelet power 
spectrum, a^, of the stars of our data set as a 
function of (cr(%)), i.e. the mean value of their 
cr-spectrum. Stars for which no cr-spectrum could 
be computed have been removed. For each star, 
marked in the plot by its WR number, a la error 
bar is plotted in both axis. The vertical dotted line 
marks the limit of cr > 5% adopted in this study, 
and the vertical one marks > 250 km s~^. All 
LSV stars are in the upper-right quadrant, except 
for WR63 (see text). 



8 LSV. Among WNL stars (including WN6 stars 
with hydrogen), 2 are NV, 10 SSV and 8 LSV. In- 
terestingly, all the WN7 stars without hydrogen 
present in our sample are, with the WN8 stars, 
the only WNL stars classified as LSV. Due to the 
proximity of the two spectral types, it remains to 
be seen if the variability observed in WN7 stars 
is similar to the variability observed in WN8 stars 
and possibly related to pulsations. Only a bet- 
ter monitoring of these stars with a more inten- 
sive temporal sampling than that obtained for this 
study can provide more clues towards answering 
this question. 

Concerning spectral variability, we find no clear 
difference in the fraction of stars we find in the 
diffcrcint categories among WNE and WNL stars. 
Indeed, we find 18% 45% and 36% of NV, SSV and 
LSV respectively amoung WNE stars compared to 
10%, 50% and 40% of NV, SSV, LSV respectively 
amoung WNL stars. For stars showing small-scale 
variability, the numbers above may be mislead- 
ing. Indeed, small-scale variations could be eas- 
ily masked by large-scale changes. Therefore, it 
is impossible from the present study to provide 
an accurate fraction of stars showing this tyj)e of 
variability and the numbers above represent, in 
that case, only lower limits. For stars present- 
ing large-scale variations, the difference between 
WNE and WNL stars is not in the fraction of stars 
we classified in the different categories, but rather 
in their distribution among the subtypes. Indeed, 
while among WNE stars almost all spectral sub- 
types show some stars with large-scale variability, 
within WNL stars there is a striking dichotomy: 
None of the 9 WN6h, WNTh stars show large-scale 
variability while all 8 WN7 and WN8 stars do. 

Among WC stars, only those of spectral type 
WC9 show LSV. The only WC9 star which is not 
a LSV is WR103, which we classified as SSV and 
for which the cr-spectrum barely reaches the 5 % 
level, i.e. just at the limit to be classified as LSV 
according to the criteria we adopted in Paper I. It 
has been proposed in Paper I that the high level 
of variability of WC9 stars could potentially be 
related to excess emission from a yet unidentified 
shock cone formed following the collision of two 
stellar winds in a binary system. Note that in prin- 
ciple, the stars in our sample are not known to be 
binaries and this remains to be checked via a more 
intensive spectroscopic monitoring. Interestingly, 
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all the WC9 stars in our sample show a high level 
of variability, even if they are not known to be dust 
producers, such as WR81, WR88 and WR92. 
Moreover, the only dust-producer of WC8d sub- 
type in our sample, WR53, is classified as SSV, 
since its a-spectrum reaches only a level of 2 and 
3 %. Hence, we conclude that the large-scale vari- 
ability is a common feature of all WC9 stars, but 
not necessarily of all dust-producers. What the 
origin of this variability is, either CIR-type or ex- 
cess emission from a wind-wind collision region, 
remains to be identified. This can only be done 
through a more intensive spectroscopic observing 
campaign; only then can the kinematical signa- 
ture of binary-related or CIR-type changes be rec- 
ognized. Ah other WC stars are NV (7 WCE) or 
SSV (3 WCE and 7 WCL). None of the WCE stars 
in our sample shows large-scale spectral variabil- 
ity. Therefore, the determination of rotation peri- 
ods from CIRs seems to be impossible for this type 
of star. The reason for this remains to be deter- 
mined. However, this is unfortunate, since WCE 
stars are thought to be possible progenitors of 
the long-soft Gamma-Ray Bur sts ( Wooslev 1993L 
MacFadven fc Wooslev I Il999t iMacFadven et al 



20011 ). Even if CIRs were present in all WR winds, 
it is possible that they are not visible in all emis- 
sion lines. For example, they may remain confined 
to regions close to the surface of the star in which 
case they may not be present in line-formation re- 
gions very far from the base of the wind, which 
is the case for the flat-topped CiiiA5696 line. It 
is clear from the different line shapes of ClllA5696 
and ClvAA5802/5812 that the former is formed 
further out in the wind than the latter. However, 
it may still be too far fro the CIRs to be deteted 
in this line. 

Finally, two WN/WC stars, WR8 (WN7/WCE) 
and WR58 (WN4/WCE), are identified as LSV. 
In addition to its line-profile variability, WR58 
shows a variable dilution in the spectral lines most 
likely due to variability in the continuum. More 
observations are needed to determine if these vari- 
ations are caused by the orbit of a companion star, 
but we prefer to keep WR58 in our list of CIR- 
type variable candidates for now. The only other 
WN/WC star in our sample, WR7a (WN4/WC), 
is classified NV. 

We wish to emphasize once more that stars that 
we have classified as NV are not necessarily con- 



stant; it is simply that no variability has been de- 
tected in spectra with a given noise level. A more 
stringent limit for the absence of variability can 
be established for some stars. For example, in the 
case of WR2, we have obtained 5 spectra with a 
typical SNR — 120 from which we deduce from 
the E-spectrum that the variability level is lower 
or equal than ~ 0.8 % of the HeiiA4686 line inten- 
sity. Moreover, we cannot detect changes in WR2 
if it has a timescale a lot shorter than the expo- 
sure time (~ 30 mins) or greater than the time 
spanned by our time-series (96 days). This estab- 
lishes a rough limit of detectability of variability 
in the case of WR2. Of course, due to our poor 
sampling, it is clearly not enough to cover all the 
possible periods/timescales in the time range span 
and all the 5 spectra could have been taken in a 
"quiet" phase of the star. In Table [21 we present 
upper limits for the spectroscopic variability for 
all NV stars in our sample, except for WR21a 
which does not show emission lines at all in the 
chosen spectral range. We list the star's name, 
the emission line chosen to carry out the analysis, 
the maximum relative amplitude of the variability 
compared to the line flux, and the total exposure 
time range spanned by the observations in each 
case. The line chosen here is the strongest, non- 
blended line available in the spectral range. 



Table 2 


: Upper limiti- 


for the level of variability 


in 


NV stars 








WR spectral line 


max. rel. 


exp. time time range 






amplitude 


(min) 


(days) 


2 


HeiiA4686 


0.8 % 


30 


96 


4 


CivA5016 


1.2 % 


30 


120 


5 


CivA5016 


1.2 % 


30 


92 


10 


HeiiA5412 


1.5 % 


30 


5 


15 


CiiiA5696 


1.2 % 


30 


9 


18 


HeiiA5412 


1.2 % 


30 


40 


33 


CiiiA5696 


1.3 % 


30 


4 


52 


CiiiA5696 


1.5 % 


30 


5 


111 


CiiiA5696 


1.1 % 


10 


30 


154 


CivA5016 


1.2 % 


30 


483 



From all LSV stars identified, we have estab- 
lished a list of promising new CIR-type variable 
candidates. Of course, from this list, we removed 
all the WN8 and WC9 stars, since it is possible 
that they are highly variable for other reasons than 
the presence of CIRs in their wind. Therefore, the 



12 



10 remaining candidates are WRl (WN4), WR8 
(WN7/WCE), WR44 (WN4), WR55 (WN7), 
WR58 (WN4/WCE), WR61 (WN5), WR63 
(WN7), WRIOO (WN7), WR115 (WN6) and 
WR 120 (WN7). Including the two already known 
cases, WR6 (WN4) and WR134 (WN6), we get 
12 CIR-type variable stars out of a total of 68 
sampled WR stars (i.e. ~ 18%). If we exclude the 
10 WN8 and WC9 stars from the total number 
of stars in the sample, this fraction of CIR-type 
variables becomes 12/58 ~ 21%, a fraction which 
is equivalent to the fraction already determined 
in Paper I, but now based on a more complete 
sample. Note that the WN4h star WR46 also 
shows clear, rapid, large-scale variability in pho- 
tometry and spectroscopy. X-rays observations do 
not contradict the possibility of the presence of 

JImil), 



□ ssv 



□ NV 



CIRs (Gosset et al. 



however intense mon- 



itoring campaigns ha ve l ed IVeen et al. I (120021) 



Oliveira et al. I ( 2004 ) and iHenault-Brunet et al 



(|201l[ ) to conclude that the origin of the changes 
is most likely non-radial pulsations. Indeed, for 
that star, the period determined in the various 
datasets, although always of the same order of 
magnitude (~8 hours), is not always strickly the 
same. This is highly incompatible with a rotation 
origin of the variability. 

Using the stellar and atmospheric parameters 
of our 10 n ew CIR-type variable candidates de- 
termined bv lHamann et all (|200d ) from a quanti- 
tative spectral analysis using their Potsdam WR 
blankeded atmosphere model, we attempted to de- 
termine if these stars present any particular physi- 
cal characteristic that could reveal information on 
the nature or origin of this type of structure in 
the wind. To do so, we have compared the dis- 
tribution his tograms of our 10 can didates to that 
of the entire iHamann et al.l (|2006l ) sample for the 
following parameters : mass-loss rate, wind termi- 
nal velocity, luminosity, radius and temperature. 
We performed a cross-correlation between the two 
distributions to check for differences and nothing 
striking was found, although our results are based 
on only 10 candidates. Ideally, this should be veri- 
fied once a larger CIR candidates sample has been 
identified. 

The next step is to determine whether our can- 
didates really show CIRs. To do this, we have 
initiated more intensive monitoring of our targets 
in photometry and in spectroscopy on different 



Fig. 24. — Histogram of the spectral type of the 68 
studied stars. This includes the 25 stars discussed 
in Paper I, the 39 stars presented in this paper and 
4 stars for which spectra has alread y been pub- 
lished in the lit erature, i.e . WR6 (IMorel et al 



19971). WR40 (lAntokhin et al. I Il995l). WR46 



(IVeen et al. 1 |2002|) and WR 138 (fLepine et al 



19961 ). Marked in black are the stars identified 
as LSV, in grey the stars we classified as SSV, and 
the NV are in white. 



timescales ranging from 10 mins to several days 
over a long period of time in order to search for 
periods in the large-scale spectral variations. To 
date, we have obtained several nights on different 



telescopes to monitor WRl (jChene fc St-Louis 



2010l ). WR115 and WR120. Also, we must keep 
in mind that, even if we classify them as SSV, 
WR71 (WN6), WR75 (WN6) and WRl 10 (WN5- 
6) also show interesting cr-spectra, and they have 
to be considered as second priority targets. Fi- 
nally, we believe it would be extremely useful to 
extend the sample of observed stars to fainter tar- 
gets in order to identify more CIR candidates and 
thus generate a more statistically viable sample. 
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